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Abstract. In this paper, we investigate the characteristic structure of 
the full equations of magnetohydrodynamics (MHD) and show that it 
satisfies the hypotheses of a general variable-multiplicity stability frame- 
work introduced by Metivier and Zumbrun, thereby extending to the 
general case various results obtained by Metivier and Zumbrun for the 
isentropic equations of MHD. 



1. Introduction 

Stability of hyperbolic boundary value problems and the related shock 
problems has been studied by a number of authors for years. There are 
satisfactory theories for symmetric systems with dissipative boundary con- 
ditions or hyperbolic systems with constant multiplicities provided that the 
boundary conditions satisfy Lopatinki conditions. Many applications fall 
into one of these two cases. But the equations of magnetohydrodynamics 
(MHD) fail to satisfy the assumption of constant multiplicities. Moreover, 
though they are symmetrizable, boundary conditions arising in the shock 
stability problem are not dissipative. This motivates the further study of 
the variable multiplicities case. It is proved by [MeZ] for symmetrizable 
systems that under certain conditions on eigenvalues, the uniform Lopatin- 
ski condition is equivalent to the maximal energy estimate which gives lin- 
earized and nonlinear stability. Specifically, they extend the construction 
of Kreiss' symmetrizers to the case that eigenvalues are not geometrically 
regular provided that they are "totally nonglancing" or "linearly splitting" 
in a sense that they define. (See section 3 for their definitions) This yields 
the maximal estimate when the uniform Lopatinski condition is satisfied. In 
[MeZj . it is also showed that all characteristics of isentropic MHD are ei- 
ther "geometrically regular" or "totally nonglancing" . Hence, a satisfactory 
framework for the study of linearized and nonlinear stability of isentropic 
MHD is achieved, through which various results are obtained. The goal of 
this paper is to investigate the eigenvalues of full MHD equations. It turns 
out that all characteristics are algebraically regular. More specifically, they 
are either geometrically regular or totally nonglancing as the isentropic case. 
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As a consequence, all results obtained in [MeZ] for the isentropic case extend 
immediately to the full MHD equations. 

2. Definitions 
Consider a first-order system 

d 

(2.1) L(p, dt, d x )U = U t + J2 A*{p)U Xj , 

i=i 

and its symbolic operator 

d 

(2.2) L(p,r,0=rld + J2^(p). 

i=i 

The characteristic polynomial of (|2.2p is P(p,r,£) := detL(p, r, £). 

Definition 2.1. (PI) T/ie system (12. ip is hyperbolic if X^=i £jAj(p) has 
all real, semisimple eigenvalues for all £ £ 

(D2) T/ie system (|2.1I) is symmetric hyperbolic in the sense of Friedrichs 
if there is a positive definite smooth symmetrizer Sip) such that S(p)Aj(p) 
is symmetric for all j = 1, ...,d. 

Let (pO)T"0)^o) be a root of P(p,r,£) = whose multiplicity is m in r. 

Definition 2.2. (a) (po> r o>£o) *s algebraically regular if on a neighborhood 
ijj of (po,£,o), there are m smooth real functions \j(p,£), analytic in £, such 
that for (p,£), 

(2.3) P(p,r,0 = e(p ) T,e)n^i(r + A i (p,0) 

where e is a polynomial in r and smooth coefficients in £ and p suc/i £/iai 
e(Po,r ,£o) / 0. 

(6) (poj^Oj^o) «s geometrically regular if in addition there are m smooth 
vectors Vj(p,£), analytic in £, swc/i i/iai 

(2.4) A(p,0vj(p,0 = ^(P,0vj(P,0 
where vi,...v m are linearly independent. 

Definition 2.3. A root (po> r Cb %> Co) of P, of multiplicity m inr is nonglanc- 
ing if the m th order of Taylor expansion of P at (po,tq,t]q,^q) satisfies 

(2.5) P m {dx) + 

where dx is the (space-time) conormal to the boundary. Moreover, it is called 
totally nonglancing if it is either totally incoming or outgoing. 

Assumption 2.4. Suppose that L is symmetric hyperbolic and nonchara- 
teristic with respect to the boundary xj = 0, and that its characteristic roots 
are geometrically regular or totally nonglancing. 
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Consider a hyperbolic system 

(2.6) L(p,d t ,dz,d Xd ) = A d (p)(d Xd + G(p,d t ,dz)) 

with the assumption that the boundary {x d = 0} is noncharacteristic, i.e., 
det A d (p) + 0. 

The classical plane wave analysis yields the boundary value problems 



(2.7) 
where 



L(p,()U = d Xd U + iG(p,()U = f, 



d-l 

(2.8) G(p,C) = A-\p)([t - i 7 )Id + J2vjA j (p)) 

i=i 

with a new variable C = ( T — *7> v)- 

Let E_(p, () be the invariant subspace of G(p, () associated to the eigen- 
values in {Im (j, < 0}. 

Definition 2.5. The Lopatinski determinant associated with {L,M){p,Q) is 
defined by 

(2.9) D(p,C) :=det(2?_(p,C),kerM(p,0). 

We say that (L, M)(p,Q satisfies the uniform Lopatinski condition on a 
neighborhood uj if there is a constant c > such that 

(2.10) \D(jp,0\>cforV( P ,OecJxSi. 

3. Computation 
The inviscid full MHD Equations are 
p t + div(pu) = 0, 

p(u t + u ■ Vu) + VP- curLB xB = 0, 
pe e (9 t + u-V6) + 6P e divu = 0, 
{ B t + curl(£ x u) = 
with the constraint 

(3.2) divP = 0, 

where p S R is density, u G R 3 is velocity, -P(p, 6*) S R is pressure, BeM 3 
is magnetic field and 9 £ R is temperature. 

Using the identity curl(P» x u) = {d\vu)B + u ■ V-B — B ■ Vu — (divB)u and 
the constraint (|3.2p . (|3.ip can be put in the symmetric form 



(3.1) 



-1 
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pt + u- Vp + pdivu = 0, 

ut + u ■ Vu + p" 1 (P p Vp + PgVfl) - p- 1 (curlS x B) = 0, 

(3-3) < 9Pf> 

6 t + u-Vd+ — divu = 0, 

pee 

B t + u-VB + (divu)£ - B ■ Vu = 0. 
Linearizing (|3.3j) about (p,u,9, B), we have the linearized equations 
Ap + pdivu = 0, 

V 



(3.4) 



A« + p~ l (P p Vp + P e V6) - p- 1 (cmlB x B) = 0, 

A0 + —divu = 0, 

peg 



, AB + (divu)^ — B ■ Vu = 

where Dt = dt + u ■ V is the material derivative. 

By taking Fourier Transform in x variable, we have the symbol matrix 
associated to (|3.4p 

(3.5) 



i 

/ p£ \ 

P" 1 ^ p- 1 ^^ p-\e r B-(B-(i)I) 

p- l P e e/e e i 

\ B tr t-{B-£)I / 

So, we can express the linear part of (|3,4p in the matrix form 

(3.6) 

( f>\ ( \ 

P^PppP + P~ l Pek tr + P~\B ■ B)e* - p- 1 ^ • OS* 

\ ■ ii) — (B ■ 0«* / 



\5/ 

Now, consider the eigenvalue equation associated to ([37 

/ ? \ / \ 

(3.7) A(E7,£) , =A 



\5/ 



V 5 / 



Equivalently, we have 
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(3.8) 



p(£ • u) = Xp, 

P^PpP? + P^Pek* + P~\B ■ m l - p-\B ■ £)B l = An 
p^PgO/egit ■ u) = xe, 
k (£ • u)B* - (B ■ £)n* = XB. 



We adopt the variable £ for the spatial frequencies and let 
(3.9) £ = |£|l, u\\=i-u, u± = u-u\\i. 

Then we can rewrite the eigenvalue equations f|3.8|> as 



(3.10) 



' pu\\ = V, 
P pP + P e 6 + B ± -B ± = Apu|| 
-Hq 1 B\\B± = Xpii±, 
P e 9/e e u\\ = pX6, 
= Afl||, 
[ uuB ± -Buu ± = XB ± . 



Let a := p/p and u := Bjyfp. The characteristic system (I3.10P reads 
ttn = Ad, 

Ppb + p P$0 + v± ■ v± = A?i|| , 

—v\\v± = Xu±, 

(3-11) { PeO. ~. 

nil = X6, 

e-eP 

= Any, 

nj_nii — v \tu± = Xv±. 
Now, we introduce new variables 



(3.12) 



x := i> + p- x P B e\ 
y := P e 6a - egp0. 



Note that the transformation (<j, I 



(x, y) is invertible if and only if 



pP p eg + p 7^ 0. So, it is invertible since P p + > 0. Using new 
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variables (x,y), we have the equivalent system 



(3.13) 



/re 



Ax, 



= Ay, 

x + v j_ • v± = Xuu , 
-v\\vx = An ± , 



= \v\\, 

V±U\\ — V\\U± 



Note that the zero-eigenspace of A is Eo := {x = 0, u = 0, v± = 0}, which 
is of dimension 2. On the space Eq, ^4 = 0. This implies A = u ■ £ on Eo- 
Thus, it • ^ is an eigenvalue of A of constant multiplicity 2. By inspection, 
it is both algebraically and geometrically regular. 

Taking a basis of such that v± = (b, 0) and letting a = v\\, we can 



(3.14) 



Aq + XI 



es of (i 


' ^11 




>«i 






y,v 


)' 


/A 


-7 2 

















()\ 


-1 


A 








-6 
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A 





\o 




















V 



where 7 2 



P 2 £3 

^ + > 0. 



So, the associated characteristic polynomial P{x) is 



(3.15) 
Thus, 



P(x) 



2/ 2 

X (x 

2/ 2 

x (x 



a 2 )((x 2 



a 2 )(x 2 
a 2 )(x 2 



(P, + §^))-6V) 



7 



6 2 x 2 ). 







' A 2 


= 0, 






A 2 


= a 2 


(3.16) 


< 


A 2 


= 4 






A 2 


= c 2 



/i 2 ) + V(7 2 " ^ 2 ) 2 + 46 2 7 2 ) , 
k(7 2 + ^ 2 ) " vV " /i 2 ) 2 + 46 2 7 2 ) 



where h 2 = a 2 + b 2 = \B\ 2 / p. 
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Remark 3.1. Coordinate change (13.120 is the key step in the computation, 
leading to block decomposition (13.140 , 

4. Results 

In the analysis of eigenvalues of (|3.ip . what is left is the same as the 
corresponding steps carried out in [MeZj for the isentropic case. Specifically, 
in the isentropic case treated in [MeZj, the matrix corresponding to f|3. 14j) 
has the same form, with identical 6x6 upper lefthand block and lower 
lefthand block of the same form XI^, but with dimension k = 1 instead of 
k = 2 as in the present case. Thus, decomposing into these two blocks, all 
the analysis goes as before. We summarize for completeness. 

Lemma 4.1. Let j 2 = P p + -Jr^ > 0. The 8 eigenvalues are 

Ao = £ " u (°f multiplicity 2 ), 

A±i = £-it±c s |£|, 

A± 2 = e • « ± (£ • B)/y/p, 

. A ±3 =£-U±Cf\£\ 

^((7 2 + h 2 ) + ^( 7 2_/ l 2 ) 2 + 46 2 7 2 )i 

i(( 7 2 + h 2 ) - V( 7 2 -/ i 2 ) 2 + 46V), 

/i 2 = l£l 2 /p, 

6 2 = ll x B\ 2 /p. 

Lemma 4.2. Assume that < |5| 2 / p7 2 w/iere 7 2 = P p + > 0. 
faj When £ ■ B ^ and ^ x B / 0, t/ie eigenvalues A±i,A±2 and A±3 are 
simple and Ao «s semisimple of constant multiplicity 2 (geometrically regu- 
lar). 

(b) On the manifold £ ■ B = 0, £ ^ 0, £/ie eigenvalues A±3 are simple and the 
multiple eigenvalues Ao = A±i = A±2 ar ^ geometrically regular. 

(c) On the manifold £ x B = 0, £ ^ 0, i/ie multiple eigenvalue Ao «s semisim- 
ple. When \B\ 2 < p 7 2 (resp. \B\ 2 > /ry 2 ), A±3 (resp. A-tiJ are simple and 
A + 2 = A_|_i (resp. A+3) 7^ A_2 = A_i (resp. X-3) are double and algebraically 
regular, and not geometrically regular (totally nonglancing, provided that 
u 3 -a^±B 3 /^p). 

Proof. The characteristic polynomial is the same as the one for the isentropic 
case, except additional factor x. The zero-eigenspace is of dimension 2, 
independent of £. Thus, the eigenvalue Ao = £ • u is geometrically regular. 
For other eigenvalues, the result follows from the calculation in jMcZ]. □ 



(4.1) 

where 

(4.2) c f = 

(4.3) c s = 

(4.4) 

and 

(4.5) 
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Corollary 4.3. For the linearized equations of (full) MHD, the uniform 
Lopatinski condition is equivalent to the uniform stability estimate (implying 
linearized and nonlinear stability). 

Corollary 4.4. Under the generically satisfied conditions < l-B^I 2 7^ 
p ± (7 ± ) 2 ; the uniform Lopatinski condition implies linearized and nonlinear 
stability of noncharacteristic Lax-type MHD shocks. 

As an application, we recover for the full equations of MHD the result of 
Blokhin and Trakhinin that MHD shocks are stable in the small-magnetic 
field limit B -> 0. 

Corollary 4.5. In the B — > limit, Lax-type MHD shocks approaching a 
noncharacteristic limiting fluid- dynamical shock satisfy the Lopatinski condi- 
tion if and only if it is satisfied by the limiting fluid- dynamical shock (i.e., the 
shock satisfies the physical stability conditions of Erpenbeck-Majda |Erl Maj | , 



in which case the MHD shocks are linearly and nonlinearly stable. In partic- 
ular, for an ideal gas equation of state, they are always stable in the B — > 
limit. 

Remark 4.6. Verification of the uniform Lopatinski condition has been 
carried out numerically or analytically for several interesting cases; see [BT1 . 
IBT2] . 
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